Fermions on the brane in 6D with nonsingular exponential scale factors 
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We introduce new realistic brane-solutions with exponential scale factors in the 6D-space-time. 
We show that for these solutions the zero modes of all bulk fields are sharply localized at different 
positions on the brane and have "Gaussian shape" wave- functions in the extra space. We also 
explicitly show that in the model there are cases when exactly three fermion generations naturally 
arise only through gravity. Because of localized fermion modes are also stuck at different positions 
in the extra space, there is possibility to provide a framework for natural explaining the fermion 
mass hierarchy in terms of higher dimensional geography. 
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>. . I. INTRODUCTION 

o . 

The idea of extra dimensions [1, 2, 3, 4, 5, 6] is one of the most attractive ideas concerning unification of gauge 
, fields with general relativity and new solutions to old problems (smallness of cosmological constant, the origin of the 
' hierarchy problem, the nature of fiavor, etc.). In theories with extra dimensions our world is associated with a brane, 
\ embedded in a higher-dimensional space-time with non-compact extra dimensions and non-factorizable geometry, 
i The key ingredient for realizing the brane world idea is localization of various bulk fields on a brane by a natural 
^ ' mechanism. The gravity is known to be the unique interaction having universal coupling with all matter fields, so it is 
important to find a purely gravitational trapping mechanism. The brane solutions and matter localization mechanisms 
^ I ' has been widely investigated in scientific literature [7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24]. 
, Recently the gravitational trapping of the zero modes of all bulk fields was realized for the brane solutions with 
' non-exponential warp factors in 6-dimensional bulk space-time [25, 26, 27]. In our previous article [28] in the case of 
\ the non-exponential scale factors we also have suggested new purely gravitational mechanism explaining the origin of 
■ three generations of the Standard Model fermions and explicitly show that localized fermions are stuck at different 
' points on the brane in the extra space. In these models, with non-exponential scale factors, the zero-mode solutions 
. of the bulk fields are normalizable and integrals of there Lagrangians over the extra coordinates are finite, but their 
O ^ wave functions spread rather widely in the bulk owing to the lack of the non-exponential warp factor. Thus, in order 
D , not to contradict the strict experiments such as the charge conservation law, some parameters in this models must 
' be chosen in a proper way. In some cases it is hard to say whether there is such a suitable choice of the parameters. 
^ I From this point of view it is interesting to consider the models with exponential warp factors when the wave functions 
of localized fields are sharply peaked on the brane in the extra space. At present it is known (see for example [12, 17]) 
that in 51? and QD models with exponential warp factors it is possible to localize all spin fields if one introduces 
^ : non-gravitational interactions. So it is of great interest to realize a purely gravitational localization mechanism in the 
■ - ' models with exponential scale factors. 

In this article in the case of (1 -I- 5)-space-time we introduce new brane-solutions with exponential scale factors 
sharply localizing the zero modes of all local fields (spin scalar field, spin I/2 spinor field, spin 1 gauge field, spin 
"^12 gravitino field, spin 2 gravitational field and totally antisymmetric tensor fields) on the brane in the extra space, 
and for the realistic physical situation of matter distribution in the 6D-space-time we realize the purely gravitational 
trapping mechanism explaining the origin of three generations of the Standard Model fermions from one generation 
in a higher-dimensional theory and explicitly show that fermions are sharply stuck at different points on the brane in 
the extra space. 
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II. MAIN NOTATIONS AND EQUATIONS 



The action and Einstein's equations which we consider in this article have the form: 

Rab - ^9abR = {^gAB + Tab) , (1) 

where Af , R, A, L, Rab and Tab are respectively the fundamental scale, the scalar curvature, the cosmological 
constant, the Lagrangian of matter fields, the Ricci and the energy-momentum tensors. In the bulk space-time we use 
the following metric ansatz [24, 25, 26, 27, 28] : ds^ = (jy^ {r) gapdx" dx^ — g (r) [dr^ + r^dO"^). The source of the brane 
is described by a stress-energy tensor Tab, its nonzero components we choose in the form T^fj ~ —gapFo (r) , Ty = 
—gijF (r), where Fq and F are source functions, which depend only on the radial coordinate r in the extra space. In 
this case the Einstein's equations take the following form 

(/) (j)^ ret) 2g 2^2 + 2rg (j)^Mj, ' ^ ' 



S 



d^xy/^ 



(R + 2A) + L 



Qcj)'^ 2q^'g' icf>' ^ g{F^A) 2gKpHys H" 2c^'g' _ g {F - A) 2gAphys , . 

02 ^ <f)g ^ r4> M4 ^ (/)2M2 ' c/) 02 (P^Ml, ' ^ ' 

where the prime denotes differentiation d/ dr. The constant Aphys represents the physical four-dimensional cosmolog- 
ical constant, where R^^^ — ^gapR'"^^ = -fj¥-9ai3- In this equation i?^"'^, i?*-^-* and Mp are four-dimensional physical 
quantities: Ricci tensor, scalar curvature and Planck scale. 



III. NONSINGULAR BRANE SOLUTION WITH EXPONENTIAL SCALE FACTORS 

In the case A^hys = from the equations (2)-(3) we can find [22, 23, 28] : F' 4^ (F - Fq) = ; g = 

^ ; ^^+3r^4-^== 217^ (F — A), where 5 denotes the integration constant with the units of squared length. 
The possible two solutions to this equations arc (in what follows we call these solutions "internal solutions" and 
introduce subscript "int"): 

1. The first solution with scale factor (jji^t monotonically decreasing from the value equal to 1 at the origin 2 = 
in the extra space and asymptotically approaching the the nonzero positive value equal to 1 — d at the radial 
infinity z = +oo 



1-d I 



(z) ^ Mpe-^z^-^e-^ , < d < 1, p > 2, (5 > 0, (4) 



2. The second solution with scale factor (/jjnt monotonically increasing from the value equal to 1 at the origin z ~ 
in the extra space and asymptotically approaching the positive value equal to 1-1- c? at the radial infinity z ~ +oo 

(r) = Sdpe-^zP-^e-""" , d> 0, p>2, 6 > 0. (5) 



0int (z) = 1 + d 



In (4)-(5) we have introduced the dimcnsionless radial coordinate z = r/s, where s is some positive constant with the 

units of length. In both cases the scale factor (f>i^t has an inficction point rinn.tp — £ (^F^) (^^^^^ value characterizes 
the order of the "effective width" ABianc of the "thick brane" in the extra space). These solutions have equal scale 
functions g-mt.- In the case p = 2 this function monotonically decreases from value equal to at the origin z = to 
the value equal to zero at the infinity z = +oo in the extra space and its derivative g[^^ is equal to zero at the origin 
and at the infinity in the extra space. In the case p > 2 this function is equal to zero at the origin r = in the extra 
space and monotonically increases in the range < r < rmax.g, approaches the maximum value at the r^ax.g and 
then monotonically decreases in the range Tmax.g < r < +oo and approaches the asymptotic value equal to zero at the 

radial infinity, where rmax.g = £ (^^) ' niust be mentioned that the scale factor gintif) defines the topology of 
the brane. These solutions exist for any A and corresponding source functions Fq and F can be easily calculated from 
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(40-3c)ln + 



Fo = F 



> dF 



corresponding Einstein's equations using (4) and (5): F ~ K 
where c = 1 — d. 

Similarly, as it was done in our article [2n] for the case of non-exponential scale factors, it is easy to explicitly show 
that these solutions localize the zero modes of the following bulk local fields: spin scalar field, spin 1 gauge field, 
spin '^/2 gravitino field, spin 2 gravitational field and totally antisymmetric tensor fields. The zero modes of these 
fields will be localized at different positions on the brane in the extra space, but as against a former case [ll8], in this 
case we shall have sharp localization of the modes with " Gaussian shapes" of wave- functions in the extra space due to 
exponential type of scale factors. Here we want to consider in more details the localization of the bulk spin-l/2 fermion 
in the background solution (4) (for the solution (5) the results are the same). Of course, in due analysis, in what 
follows we will neglect the back-reaction on the geometry induced by the existence of the bulk fields, and, without 
loss of generality, we will take a flat metric on the brane. Fermion action and corresponding equation of motion are 



M , 



(6) 



Introducing the vielbein through the usual definition gAB = ^^a^bVab '^ti^i'^ A,B,... denote the local Lorentz 
indices, it is easy to find the non- vanishing components of the spin-connection for the background (4): lu^" = 

— ■ Therefore the covariant derivatives have form Da'^ ~ {pa + \<^^^lzla) , = dz'^, De'^ = 

i^de + ^(^g^lzle^ ■ The 6-dimensional spinor '5 (x*^) can be decomposed on a 4-component spinor ip (x^) and 2- 

component spinor £^{z,9) : 5* (x^^) — ip (x^) (l~)£,{z,9). We require that the 4-dimensional spinor satisfies the 
equation ^^d^-il) [x^^ = mtjj (x'^) . The 2-component spinor can can be expanded over the states with fixed rotation 
momentum I in the 9 direction: ^{z,9) =J2^i (z) e^^^ , where (z) — {vi,Q) is 2-component spinor, vi and (i are 



functions of variable z. Taking gamma matrices 7^ and 7* in the form 7^ = 
we obtain the following equations for vi and Ci 



I 

1 



7 = 



-i 

1 



, as a result 



'-vi 



me /- _i_ 



/9int 



a. + 2|^ 



= 0, 

= . 



(7) 



In the zero-mass mode case m = (the case m 7^ we will consider later on) the solution to these equa- 



tions reads vi (z) = ai^^Jg-^^^' z 2+' , Q{z) = bi(j).Jg.Jz 2 ' , (z) = (jjjg.^jz- 
bi being the integration constants. Using explicit form (4) of the scale functions 0i,it 



biz' 
and gin 



, with ai and 
the normaliza- 



1 T'^ ^ ^ 

tion requirement becomes 1 = 2n {5dpe^) '^ J [afz'^'' +bfz~'^'''^ z^^^^e^^^ dz, and it is easy to see that the extra 



part of the 6I?-wavefunction is normalizablc (i.e., the integral over z converges) only in the following cases: I) 
ai ^ , k ^ , I > -p/4 ; U) ai ^ , bi ^ , I < p/4 ; III) a, ^ , 5, ^ , ~ p/i < I < p/4. 
Finally for the effective wave function of the zero-mass mode of 6-dimensional fermion in flat space we have 



'i/ ^ ip {x'^) (E) {Sdpe^) 



zif 2e i 



aiz^e'^^ 
biz-^e'"^ 



+ E 



iW 



In this case because of exponen- 



tial decreasing warp factors normalizability of wave function means sharp localization of corresponding mode on 
the brane. These components of zero-mass mode of 6-dimensional bulk fermion are localized at different points 

?'max.,i — £ { ^^^p~^ ^ ' ill the extra space. Each component can be regarded as 4-dimensional massless fermion. So in 
this case there are infinite number of massless fermions on the brane. This situation indicates that in order to single 
out finite number of these localized fermion components we need some new mechanism, and the one possible way we 
will consider in the next section. 

To conclude this section let us consider the solution to Einstein's equations without any sources {Fq (r) = F (r) = 0) 

in the case A 7^ 0, h-phys = 0. In this case the Einstein's equations reduce to -I- 3r^ If ^ ~2mt^ ' 9 = 
— , 6 — const. The solution to these equations is (in what follows we call it "external solution" and introduce 
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corresponding subscript "cxt") 

, , 1 , . Sb 1 , SA 

(Pe^t [Z) ^ —J- , 5oxt ^ = ^, 0=——, (8) 

where the a is some new positive parameter and the constants e, S and the radial coordinate z are the same as above. 
This solution is nonsingular in the range < z < +oo in the extra space. Its scale factors ^ext. and ^oxt. are 
monotonically decreasing functions in this range. 



IV. THREE FERMION GENERATIONS ON THE BRANE 



Now imagine the following realistic physical situation. Suppose in the extra space in the range < z < zq {in 
what follows we name this area in the extra space by core) we have nonzero source functions Fq (z) and F (z) which 
correspond to the solution (4) with decreasing scale factor (/)int., and out of the core, i.e. in the range z > zq there 
are no source functions, i.e. Fq (z) = F (z) = if z > zq- At present we do not argue the origin of such matter 
distribution in the 6-D-spacc-time. Let us denote by (/) and g the scale factors of solution to Einstein's equations in 
this special case. It is obvious that they can be constructed from the scale factors of solutions (4) and (8) in these 



two regions in the following way: (j) (z) 



hnt (z) , 
icxt {z) , 



:iz) 



5int (z) , 
5oxt (z) , 



< Z < Zq, 



Z > Zq. 



To avoid 



< Z < Zq, 

z > zo, 

singularities we impose on the scale factors and there first derivatives the continuity conditions at the boundary of the 

core z = Zo : 0int (^o) = 0cxt (^o) ; 0|nt (^o) = ^oxt (^o) ; 5int {zq) = ffcxt (^o) ; 5int (^o) = 5cxt (^o) • So far there were the 
following free parameters in the model: M, A, S, d, s,p, a, zg. Solving the continuity conditions we can fix the following 



parameters a 
and A (x) = (x 



A (xo) , b 



SA 



(ao-l)e'' 



-, where we have denoted xq 



> 1 



1) + a; + 1. Now we have following free parameters in the model: M, S,e,p and zg > 1. 
Let turn our attention to the localization of the bulk spin--'^/2 fermions on the branc described by our solution. In 
this case equations for the components vi and Q of the radial part (z) of zero-mass wave-function have the same 
form as in the previous section, and the solution reads 



^i{z)^(f) Jz 2 



a;z 

hz-' 



1-2 -f -J 
"int5int ^ ' 



f'oxtSoxt Z 



aiz' 
biz-' 

aiz'' 
biz-' 



< z < zo. 



z > Zo, 



(9) 



with the following normalization requirement 



1 = 2ne 




,2 21 



afz^' + bfz 



'2l\ 



dz 



afz'^'- 



+ b^ 



'21 



z (a -|- 61nz) 



-dz 



(10) 



Analysis of convergence of the last integral at the boundaries of the whole integration domain (i.e., at the points 
z = and z = +oo) shows that normalizable components of the zero- mass mode exist only in the following cases: I) 
a; 7^ , 5; = , — p/4 < / < ; II) a; = , bi ^ , < I < p/A . So the number of normalizable zero-mass modes 
depends on the value of parameter p. Namely in the case 12 < p < 16 we have exactly three normalizable zero-mass 
solutions with I = ±1, ±2, ±3. In this case the effective wave-functions of this zero-mass modes in fiat space are: 



4' (x 



ilj{xf^)<» {Sdpe 

V'(a;^) «) ((5te2)* z-i (a + blnz) 



E 

0<i<f 

-1 



a-iz-'e-'"^ 
biz-'e''^ 



, < z < Zo, 



E 

0<i<f 



a^iz-'e-'"^ 
biz-'e'"> )'^>^o- 



(11) 



radial direction 

\i/p 



of the extra 



space at 

i/p 



the following points 



rfci-m,l = l = 



They are localized in the 

e (1 - 6/p)^^^ , rferm,z=2 = £ (1 " W/p)'"' , rform,i=3 = £ (1 " 14/p)''*' . It must be mentioned that if 12 < 6 < 14 
the last formula does not work, but in this subcase directly from the wave function (11) we have rferm,i=3 — 0. For 
the boundary of the core we have zo > 1, so each of these three modes are sharply localized inside the core and 
the wave- functions of these localized modes have the "Gaussian shapes" in the extra space. These three components 
of a higher dimensional fermion (stuck at different points on the brane in the extra space) can be identified with 
different 4-dimensional fermion generations, and so we have purely gravitational mechanism explaining the origin of 
three generations of the Standard Model fermions. 
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V. NON-ZERO-MASS MODES 



Now let us consider the case m 7^ for the physical situation described in the previous section. We have to 
solve the equations (7) in two regions - inside the core and outside the core. We will seek the solution in the 
form vi (z) = z'~2 ^ (2) , Q (z) = (j)~'^g~iz~''~2B{z), where 4> = 4'int and g = gint if < z < zg, and 



/'cxt , 9 = 5cxt if z > zq. From (7) for the functions A and B we get equations ^ 



(p ^21 dA 



m 



^ ' S "^^"^'^ — 'Ti^e^^z^^'i? = . We examine the solution A{z), the results for the B{z) are the same. Taking 
into account (8) and normalization requirement for wavefunction at the radial infinity in the extra space, outside the 
core the equation for Acxt and its nomalizable solution assume the following form: A'^^^ + + ij x^^^^ — ^^i^^^^^. = 

, Acxt [x) — ax p , where a is the integration constant and we have introduced new variable x = z^. Inside 
the core the equation for Ajnt has the form 



d 

dx 



ii+v dAint 

a [x) X 2p 



dx 



^ Aint = 0, where ct (x) = [l - d (l - e"'')] e^ . (12) 



p a {x) 



At the origin of extra space on the solution we impose condition \A\^^ (0)| < +00, and at the boundary of the core 
a; = xq we require the following matching conditions: Ajnt {xq) — A^xt (xq) , A[^^ (xq) = A'^^^ (xq). So in the core, < 
X < xq, we have a Sturm-Liouville differential equation for Ai^t with the following boundary conditions: |Aint (0)| < 
' ^^aIx^)^ ~ — '^'^pjp™'''^'' ■ In the vicinity of the extra space origin for the asymptotic form of (12) and for its 
normalizable solution we have: xA" + ^2^ + A' — bI-Ma = , Ajnt (x) — Px^ BesselJ {v,lo^/x) , where v = 

, uj — 2y^^p^ and /3 is the integration constant. To make some interesting qualitative conclusions we consider 
this solution as the correct solution on the entire core range and choose xq sufhciently large to use the asymptotic 
expression for the Besscl functions ( BesselJ{v, y)lj,^+oo ^ \/^*^°^ (2/ ~ ^ ~ f ) ) ^^"^ boundary of the core . Then 
for the eigenvalue A we have the following transcendental equation: tan ( 2A ~ y j — ^^^y ^ + ^d^^^ ' ^ ^ 



. This equation gives the discrete spectrum for A = Ai, A2, A3, ... , and for corresponding nonzero masses 



we get: m„ ~ y JdxZ -Between zero and nonzero modes there is the gap depending on the value of parameter 5, 

so choosing the appropriate value of 5 the non-zero modes can be neglected at low energies. The same conclusions 
can be done with respect to the exact solution of (12), although it is difhcult to find the solutions explicitly. 



VI. CONCLUSION 



To conclude we summarize our main results. In this article in the case of (1 + 5)-space-time we have introduced new 
brane-solutions with exponential scale factors sharply localizing the zero modes of all local fields (spin scalar field, 
spin spinor field, spin 1 gauge field, spin gravitino field, spin 2 gravitational field and totally antisymmetric 
tensor fields) on the brane in the extra space. We have explicitly shown that these solutions localize infinite number 
of 4_D-fermions on the brane. To solve this problem we have introduced another specific brane-solution with realistic 
source functions, i.e. the source functions Fq and F are nonzero inside the core < z < zg in the extra space, and 
Fq = F = outside the core. We also have explicitly shown that in this model there are cases when exactly three 
4£'-fermion generations can naturally arise only through gravity. As the zero modes of all other local bulk fields, these 
three zero modes of bulk fermion are sharply peaked at different positions on the brane and have " Gaussian shape" 
wave-functions in the extra space. So there is possibility to provide a framework for understanding the fermion mass 
hierarchy in terms of higher dimensional geography [2!)]. 
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